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■ ABSTRACT 

Winds from protogalactic starbursts and quasars can drive shocks that heat, 
. ionize, and enrich the intergalactic medium. The Sedov- Taylor solution for point-like 

explosions adequately describes these blastwaves early in their development, but as the 
time since the explosion (t — ti) approaches the age of the universe (t), cosmological 
effects begin to alter the blastwave's structure and growth rate. This paper presents 
an analytical solution for adiabatic blastwaves in an expanding universe, valid when 
J> ■ the IGM is homogeneous and contains only a small fraction of the total mass density 

(^^IGM ^o)- In a flat universe, the solution applies until the age of the universe 



5 



. approaches tiJljQj^ , at which time the self-gravity of the matter associated with the 

. shock compresses the shocked IGM into a thin shell. When Oigm ^ 0.03, blastwaves 

\^ i starting after z ~ 7 and containing more than 10^^ erg remain adiabatic to relatively 



low z, so this solution applies over a wide range of the parameter space relevant to 
galaxy formation. Using this analytical solution, we examine the role protogalactic 
explosions might play in determining the state of intergalactic gas at 2; ~ 2 — 4. Since 



^ • much of the initial energy in galaxy-scale blastwaves is lost through cosmological 

■ effects, photoionization by a protogalaxy is much more efficient than shock ionization. 

Shocking the entire IGM by 2; ~ 4, when it appears to be substantially ionized, is most 
easily done with small explosions 10^^ erg) having a high comoving number density 
^ I lMpc~^). Larger-scale explosions could also fill the entire IGM, but if they did 

I so, they would raise the mean metallicity of the IGM well above the levels observed 

in Lya clouds. Since the metal abundances of Lya clouds are small, the metals in 
these clouds were probably produced in small-scale bursts of star formation, rather 
than in large-scale explosions. The H I column densities of protogalactic blastwaves 
are much smaller than those of typical Lya clouds, but interactions between shocks 
and preexisting Lya clouds can potentially amplify the neutral column densities of 
preexisting clouds by a large factor. 

Subject headings: galaxies: intergalactic medium — galaxies: quasars: absorption lines 
— ISM: supernova remnants — shock waves 
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1. Introduction 

Gravity initiates the formation of galaxies in the early universe, but shortly thereafter, 
explosions, winds, and ionizing radiation from massive stars and active galactic nuclei complicate 
the hydrodynamics of galaxy formation and regulate the transformation of gas into stars (see 
recent reviews by Bond 1993 and Shapiro 1995). In the nearby universe we see starbursting 
galaxies expelling up to 10^^ erg of thermal and kinetic energy into intergalactic space over the 
lifetime of the burst (see Heckman, Lehnert, & Armus 1993 for a review). The explosive output 
from a protogalaxy could be even larger, if most of its stars formed on the dynamical timescale of 
the galaxy. 

Once the hot gas from a starburst escapes its host it flows into the intergalactic medium 
(IGM), shocking it, and eventually forming a blastwave similar to a supernova remnant. Because 
cosmological effects, such as the expansion of the universe, eventually govern their behavior, the 
explosions that emanate from quasars or protogalaxies have sometimes been called cosmological 
blastwaves. As these shocks pass through intergalactic space, they compress the IGM, sweeping 
it into denser sheets. Such a picture has naturally motivated suggestions that cosmological 
blastwaves are somehow related to the intergalactic Lya clouds whose absorption signatures 
striate the spectra of high-redshift quasars (Ozernoi & Chcrnomordik 1978; Ostriker & Cowie 
1981; Chcrnomordik k Ozernoi 1983; Ikeuchi & Ostriker 1986; Vishniac & Bust 1987). 

Early treatments of cosmological blastwaves posited a baryonic IGM that dominated the 
matter content of the universe (Schwarz, Ostriker, & Yahil 1975; Ostriker & Cowie 1981; Ikeuchi, 
Tomisaka, &: Ostriker 1983; Bertschinger 1983). Today the most successful cosmological models 
assume baryons comprise a relatively small fraction of the universe's matter content. Since the 
older solutions do not apply, cosmological blastwaves and their impact on the IGM deserve a 
fresh look from a more contemporary perspective. Computational approaches to such "feedback" 
during galaxy formation, while rapidly growing more sophisticated, are still limited by numerical 
resolution (e.g. Cen & Ostriker 1993; Navarro & White 1993). Analytical approaches, like the 
one taken here, while simplistic and idealized, can be more flexible than numerical solutions and 
offer insights complementary to those gained from simulations. Here we demonstrate that, in a 
universe filled primarily with collisionless matter, the structure and growth rate of a cosmological 
blastwave can be determined analytically. This new analytical solution turns out to be handy for 
reassessing the significance of protogalactic shocks in the IGM. 

The paper is organized as follows. Section 2 presents an analytical solution for adiabatic 
cosmological blastwaves, valid when the mass density of the IGM is gravitationally unimportant. 
This solution is directly analogous to the Sedov- Taylor solution for adiabatic supernova remnants. 
As long as the shock remains adiabatic, the postshock density does not exceed four times the 
ambient intergalactic density until the shell of collisionless matter that collects at the shock front 
becomes self-gravitating. Section 3 discusses how ionizing radiation modifies the medium into 
which intergalactic shocks propagate, and § 4 considers the validity of the fluid approximation 



-3- 



in the IGM. Section 5 explores the global impact of protogalactic shocks on the temperature 
and metallicity of the IGM and derives constraints on the number density and energy scale of 
blastwave sources. Section 6 examines possible relationships between intergalactic blastwaves and 
Lya clouds, and § 7 recapitulates the important results. 



2. Blastwave Development 

Intergalactic shocks can have an enormous impact on the IGM. Early calculations assuming 
that the ratio of the IGM mass density to the total mass density of the universe (Jigm = ^^igm/^^o) 

was close to unity found that protogalactic blastwaves could fill the IGM, heating and ionizing 
it by rcdshift z of a few (Schwarz et al. 1975). Subsequent treatments attempting to resolve 
cosmological shock structure in more detail (Ikcuchi et al. 1983; Bertschinger 1983; Vishniac, 
Ostriker, &; Bertschinger 1985) generally took 5igm = 1 and found that shocks in an IGM 
dominated universe formed dense shells. Ikeuchi et al. (1983) assumed that dense shells would also 
form when Sjgm ^ 1- At late times this certainly happens. An explosion in the IGM at a time 
ti since the big bang creates a cavity within which the mass density is a fraction 1 — 5iqm lower 
than the background density. In a flat universe (17 = 1), collisionless matter within the cavity 

—3/2 

gradually overtakes the shock front and forms a shell when the age of the universe is ~ ^i^igm 
(Bertschinger 1985). After this time, blastwave evolution asymptotically approaches the usual 
self-similar solutions in which ti can be neglected, and the shocked IGM collapses into a dense 
layer just behind the shock front (Ikeuchi et al. 1983; Bertschinger 1983). If 5igm *C 1, most of 
the interesting interactions of protogalactic shocks with the IGM occur long before the blastwave 
reaches the asymptotic self-similar stage. 

This section presents an analytical solution for adiabatic cosmological blastwaves valid when 

3/2 

SiGM 1 and t <^ ti^iGM • Taking advantage of some similarity properties of cosmological 
expansion (see Shandarin 1980 or Shapiro, Struck-Marcell, &: Melott 1983), one can transform 
the cosmological fluid equations into a form without an explicit Hubble flow or background 
density. When the adiabatic index, 7, of the IGM equals 5/3, the transformed fluid equations 
admit a solution analogous to the Sedov-Taylor solution for adiabatic supernova remnants in 
a homogeneous medium. This solution turns out to be widely applicable, since cosmological 
shocks at 2; ^ 7 tend to remain adiabatic. The solution breaks down when either cooling becomes 
important or perturbations in the collisionless background flow become large. 



2.1. Analytical Solution 



Self-similar solutions for cosmological blastwaves can be found straightforwardly when the age 
of the blastwave is much less than the current age of the universe or when the age of the blastwave 
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2 -i.+P'^- (2-4) 



is infinitesimally close to the current age of the universe (see Ostriker Sz McKee 1988 for a review). 

In these cases, there is only one important timescale. The intermediate case has been considered 
unfit for self-similar treatment because both the current age of the universe (t) and the age of the 
blastwave {t — ti) are important. However, when 6igm -C 1 a self-similar solution for cosmological 
blastwaves exists if 7 = 5/3 and the blastwave is adiabatic. 

To find this solution, we begin with the Newtonian fluid equations, valid as long as the length 
scales are much smaller than the horizon of the universe and the velocities are non-relativistic: 

^ + Vpn = (2-1) 

^ + (u-V)u = -V$--VP (2-2) 
ot p 

= At,G{p + pg) (2-3) 

de P dp 

di ~ ^'dt 

e = . (2-5) 

7- 1 p 

Here, r is the coordinate vector, u is the corresponding velocity, P is the gas pressure, e is the 
specific internal energy, p is the IGM mass density, and pc is the mass density of collisionless 
(non-baryonic plus condensed baryonic) matter. The undisturbed IGM comoves with the Hubble 
flow. Deviations from this flow can be expressed as a peculiar velocity field v = u — H{t)r under 
the influence of a perturbed potential fleld ^ = $ — 27rGp(t)r^/3. In terms of the scale factor 
a{t) = [1 + z{t)]~^, which obeys a/a = H{t), the mean mass density is p{t) = pQa~^{t), where 
Po = Pcr^o = SHgilQ/SnG is the mass density of the universe at the present time (to). 

A transformation originally applied by Shandarin (1980) recasts the cosmological fluid 
equations into a very useful form (see also Shapiro et al. 1983). The transformed variables are: 

dt = a''^dt , f = a~^r , v = av , p = a^p 

P = a^P , e = a'^e , ^ = a^cf) . (2-6) 
These variables obey the equations 

f + Vpv = (2-7) 

^ + (v.V)v = -V0-lvP (2-8) 
dt p 

= iTiGip + pG - Po)a (2-9) 

P 

di p'^ di ' ^~ ~ " \a di . 



- ^^^ + (5-37)f-i)6 + /5aH-p^a-^A (2-10) 



1 P 
7-1/5 



(2-11) 
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The form of these equations closely follows that of the initial set except that f and v now represent 
comoving positions and peculiar velocities. 

When (5iGM 1, perturbations in the IGM have little effect on the underlying gravitational 
potential. If we also assume that the underlying collisionless matter is distributed effectively 
homogeneously on the scale of the blastwavc, terms containing (j) can be ignored. In an adiabatic 
shock, the heating and cooling terms are negligible. The only term remaining that depends on 
cosmic time is the second term in the energy equation, which vanishes if 7 = 5/3. This term 
accounts for the change in e that arises when e does not vary as a~^. In a monatomic gas, all the 
internal energy is in translational degrees of freedom. These atomic motions can be viewed as 
peculiar velocities that redshift away like a~^, so that e oc in the absence of heating, cooling, 
and noncosmological work. Energy contained in rotational and other internal degrees of freedom 
does not redshift away and adds to e as the universe ages. Thus, in the transformed frame, the 
second term in the energy equation heats the gas when 7 > 5/3. 

The above equations are formally identical to the adiabatic fluid equations without 
cosmological and gravitational terms when 7i = A = = and 7 = 5/3. Since the time 
transformation involves only the time derivative, we are free to define t = to coincide with 
the cosmic time ti at which a galactic explosion (quasar or starburst) introduces an energy 



Eq = a '^Eq into the IGM. Solving for t{t) requires knowing a{t) and thus Qq: 



^oHq 



(2-12) 



The standard Sedov self-similar solution (Sedov 1959, 1993) then satisfies the fluid equations, 
giving a shock radius of 



Rs = 



PIGM 



1/5 



, 



(2-13) 



where ^ = 2.026, picM = Pcrf^iGM and JIjgm = ^o^iGM- At early times, i !v {1 + zi)'^{t — ti), and 
the shock radius in physical coordinates is 



Po^^igm(1 + zi] 



1/5 



(2-14) 



as expected. At large times, i asymptotically approaches a constant, so Rg approaches a constant 
comoving radius. The complete solution for the shock radius is 



32tt(,GEo (l + nozi) 



SH^QiGM f^g(l + ^i)2 



= (7.5 Mpc) 



1 + Oo-zi 
nlil + z,r 



1/5 



1/5 



1 + n^z 

1 + riozi 



pV5>-4/5p>-l/5 
-^61 "'50 ''IGM 



1/2] 2/5 



-1 



(1 + ^) 

1 + Oqz \ 



(2-15) 



2/5 



(1 + ^) 



-1 
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where Eqi = Eq / {10^^ erg) . This solution is also valid when Oigm = ^ because the 
gravitation of the collisionless matter is insignificant in this limit. 

Since JIigm 1, the asymptotic radii of cosmological shocks can exceed 10 Mpc in comoving 
coordinates. Figure 1 shows how quickly the shock radius and shocked mass approach their 
asymptotic values in a flat universe. Voit (1994), solving for cosmological blastwave evolution in 
the thin-shell approximation, obtained an identical rate of convergence to the asymptotic solution. 
At the asymptotic radius, the blastwave has spent all the input energy in lifting the shocked IGM 
out of the gravitational potential of the collisionless matter interior to the shock. In the thin-shell 
approximation, the blastwave must lift all of the shocked IGM to the maximum radius, so the 
asymptotic radius in the thin-shell case is slightly smaller than in the exact solution. 



2.2. Cooling 



The onset of efficient cooling ends the adiabatic stage of blastwave development. Inverse 

Compton scattering of hot electrons off the cosmic microwave background cools cosmological 

2/5 

shocks effectively at z > S/i^q — 1 in a flat universe. Radiative cooling is more important at 
lower redshifts but does not cool the shocked IGM very quickly. The cooling function A(T) at 
temperature T in a low-metallicity plasma is of order 10~^^ ergcm^ s~^ for 10^ K < T < 10^ K 
(e.g. Sutherland &z Dopita 1993; Schmutzler & Tscharnuter 1993). Just behind a strong shock 
front, the hydrogen number density nn is four times the ambient IGM hydrogen density and the 
gas temperature is T = SiiVg/lBk, where Vg is the shock velocity relative to the Hubble flow and 
is the mean mass per particle. When the postshock cooling time tc = 5kT/n}iA{T) is smaller than 
t, the current age of the universe, a cosmological shock is not adiabatic. If fio = 1, cosmological 
shocks cool radiatively for 



32 Ho^iGu-^ 



IStt 



Gulp 



1/2 



,3/4 



(2-16) 



A 



10 erg cm^ s ^ 



1/2 



If rio <^ z ^ , the criterion for radiative cooling becomes 



16 I^o^^igmA 



1/2 



(1 + ^) 



(2-17) 



A 



10~^^ erg cm^ s" 



\ 1/2 

T j hli' (1 + z) 



Ionization losses can become important before radiative losses do. Blastwaves propagating 
into a neutral medium expend much of their energy collisionally ionizing the neutral gas when 
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Vs 4:0 — 80 km s ^ . However, ultraviolet photons from the protogalaxy or quasar driving the 
blastwave are likely to preionize the IGM into which the shock propagates (§ 3). 

The analytical solution of § 2.1 straightforwardly gives the velocity at the leading edge of an 
adiabatic shock: 



3J^iGM (i + nozi)^/^ 



nl/5 



1 + Op-g \ 
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1 



1 + Zi 



(2-18) 



= (75 km s 



^^^(l + ^i)^ 

1 -I- UqZi 



3/10 



pl/5, 1/5^-1/5 
-^61 ''■50 ''IGM 



i + npz y/^' 
i + nozi) 
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1 + Z 
1 + Zi 



Note that Vs{z) depends only on f^o, zi, and the parameter combination E^ih^oQ^Q^. Figure 2 
shows how Vg varies with z for various values of £^61^50^igm several different initial redshifts, 
assuming = 1. The vertical line indicates where inverse Compton scattering is marginally 
effective in cooling intergalactic shocks. The effectiveness of radiative cooling depends on the 
density and mctallicity of the ambient IGM. Metallicities of Lya clouds are probably ^ 10~^ 
times solar (Tytler & Fan 1994; Cowie et al. 1995). At these levels of enrichment, A(T) is nearly 
equal to the cooling function for a metal-free plasma. Solid lines across the lower portion of 
Figure 2 show where tc ~ t, according to the cooling functions of Sutherland & Dopita (1993), for 
a metallicity of 10~^ times solar and f^icM = 0.01 and 0.03. The shocks described here do not 
cool via radiative processes as long as Vs{z) lies to the right of the Compton cooling region and 
above the radiatively cooling region. For a wide range of parameters, cosmological blastwaves are 
adiabatic. The analytical solution derived here is therefore valid in many interesting cases. 

Once cooling becomes efficient, the blastwave grows more slowly than in the adiabatic case. If 
the shock cools radiatively, but the interior does not, then Eq oc R^"^ and Rg oc p/'^ . This case is 
directly analogous to the pressure-driven snowplow solution for supernova remnants. If both the 
shock and the interior of the blastwave radiate efficiently, then blastwave momentum is constant 
in the transformed frame {Eq oc Rsi~^) and Rg oc t^/^, as in the momentum-conserving snowplow 
solution for supernova remnants. For more on supernova-remnant solutions see Ostriker h McKee 
(1988). 

Cooling reduces the ultimate radius of a cosmological blastwave and the total amount of mass 
it can shock. Let i?asy be the asymptotic comoving radius of an adiabatic blastwave, and define 
r] = i?s/i?asy Suppose that at r] = r]c and t = the blastwave begins to cool efficiently and 
obeys Rg = RasyVcii/ic)^ thereafter. The new asymptotic comoving radius of this blastwave, after 
cooling, is RasyVc ^^^"^ . Relative to the adiabatic case, the asymptotic shocked mass is r/c^'' times 
smaller when (3 = 2/7 and ifj^ times smaller when (3 = 1/4. 



2.3. Gravitation of the Shocked Shell 
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Shock compression of the IGM creates a radial perturbation in the gravitational potential 
that can modify the interior structure of a cosmological blastwavc. In transformed variables, 
the surface density of the shocked IGM shell is o" = HqQiqj^Rs/^ttG. Note that the form of 
equation (2-9) implies G = G{1 + z)~^ . Behind the shock, the gravitational pressure scale height 
in the transformed system is 

Thus, the gravitational pressure scale height in physical space is 

^9 - ;^F7^ (^) , (2-20) 



25J^iGM + 

when n = 1. In the absence of gravity, the density scale height inside an adiabatic blastwave is 
~ Rs/12, so self-gravity does not alter the radial density structure of the blastwave significantly 
until \g --^ Rs/^'2 at t ~ ti$7jQ^ . Gravitational forces modify the radial pressure profile of an 
adiabatic cosmological blastwave somewhat earlier. Without gravity, the gas pressure at r < 0.8 ii^ 
plateaus at ~ 0.4P(i?s) (Sedov 1959, 1993). The gravitational potential of the shocked shell 
steepens this plateau when Ag ~ i?s. We can always neglect the self-gravity of an adiabatic shock 
when fiiGM ^ 1. 



2.4. Fragmentation 

Density perturbations in the swept-up shell of gas behind a cosmological shock can become 
gravitationally unstable, leading to fragmentation of the shell and possibly the formation of Lya 
clouds or galaxies (e.g. Ostriker &z Cowie 1981). However, when ^igm ^ 1, gas shells behind 

adiabatic shocks do not fragment until the response of the coUisionless matter starts to govern 
the dynamics. Formation of condensed gaseous structures through shock fragmentation alone thus 
requires either shocks that cool radiatively or Siqm ~ 1. 

Let us examine the development of density perturbations in the transformed variables when 
^IGM -C = 1 and the shock is adiabatic. A density perturbation of wavelength A in a gas shell 
with a characteristic sound speed Cg grows if the sound crossing time \cj^ exceeds the free-fall 
time (A/Ga)^/2 (Ostriker & Cowie 1981; Vishniac 1983; Voit 1988). Thus, the minimum timescale, 
in transformed coordinates, for postshock IGM density perturbations to grow gravitationally is 
ig ~ Cg/Ga. While the shock remains adiabatic, Cg « 0.16 Rg/i, and ic ~ 9(1 -|- z)^}jQj^tQi~^ . 
Postshock perturbations do not grow significantly until i 3(1 -|- z)^^'^Q,^Q^tQ, which corresponds 
to the epoch in physical time when t ^ tiQ^Q^^ . At this late point in the development of the 
blastwave, we can no longer ignore the peculiar flow induced in the coUisionless matter as the IGM 
is swept up. 
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This stability analysis can be formulated more rigorously following Vishniac (1983), who 
examined the behavior of spherical blastwaves in the thin-shell approximation. Applying Vishniac's 
formalism to the gravitational stability problem in transformed coordinates yields a similar growth 
time, tg « Cs/irGa, and the same qualitative result. If (^igm 1, adiabatic cosmological shocks do 
not fragment gravitationally until the response of collisionless matter to the blastwave begins to 
drive the shock front. 

Radiative cooling, when it finally becomes a significant energy sink, reduces Cg, allowing 
gravitational instabilities to develop earlier, but isothermal shocks are subject to a disruptive 
dynamical overstability that operates much more quickly than gravitational collapse (Vishniac 
1983; Bertschinger 1986). Ripples in the shocked shell induce transverse flows of gas from the 
leading regions of the shell to the trailing regions. As the surface density of the trailing regions 
grows, the trailing regions begin to overtake the leading regions, and the process reverses. When 
the overstability becomes nonlinear, weak shocks damp the transverse flows, staving off complete 
disruption of the main shock front (Grun et al. 1991; Mac Low & Norman 1994). Gravitational 
instabilities therefore proceed to grow roughly as expected. 



2.5. Comparisons with Other Solutions 

Since the pioneering work of Schwarz et al. (1975), several investigators have presented 
numerical and analytical solutions for cosmological blastwaves. Many of these treatments assume 
^IGM = 1 and are not valid when Oigm <SC ^^o- At late times [t ^ ^i^igm ) in a flat universe filled 
with collisionless matter, 

= 1.89 (E^oG)^/^*^/^ (2-21) 

(Bertschinger 1985). If ^^igm ^ 0.1, this solution does not yet apply to blastwaves beginning after 
the era of Compton cooling (zi < 7), because the collisionless component has not yet responded 
fully to the motion of the IGM. Several workers have formulated approximate expressions for 
Rs{t) in the important regime where t <^ ^i^igm ^^'^ ^IGM <^ = 1- Here we compare these 
expressions with the exact solution derived in § 2.1. 



2.5.1. Shock Radius 

Adiabatic cosmological blastwaves start as spherical Sedov shocks and then asymptotically 
approach a constant comoving radius before the collisionless matter drives Rg to the late-time 
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solution. Ozernoy & Chernomordik (1978) constructed an approximate expression for Rg by 
matching a Sedov solution at early times to the solution of Schwarz et al. (1975) at late times. 
This approximation is poor because it presumes that the collisionless component responds to tlie 
blastwave on a timescale ~ ti. Vishniac &: Bust (1987) introduced a correction factor to the 
Sedov solution that reproduces, with reasonable accuracy, the behavior of Rs{z) in numerical 
models of adiabatic cosmological shocks (Vishniac et al. 1985) when Oigm <S = 1. Figure 3 
shows how their approximation (VB87) compares with the exact solution of § 2.1 (V95) and the 
late-time self-similar solution of Bertschinger (B85) for JIigm = 0.1 and 0.01. The corrected Sedov 
solution remains within 10% of the exact solution until t lOOti, regardless of ^^igm- After this 
time, the slope of VB87 is similar to B85. The point of departure of VB87 from V95 should be 
later for JIigm = 10~^, indicating that the VB87 correction factor was derived from models with 
^^IGM ~ 0.1. Ostriker &; McKee (1988) formulated an expression for blastwave expansion that can 
be compared more directly with the solution of § 2.1. They constructed an approximate solution 
for adiabatic blastwaves when J^igm <C J^o = 1 by matching a Sedov- Taylor solution at early times 
to a comoving shell solution at late times. The approximation of Ostriker & McKcc (1988), labeled 
OM88 in Figure 3, agrees well with the exact solution until (l-|-2:i)/(l-|-^)~4 and ultimately 
reaches a somewhat smaller comoving radius. 

The asymptotic comoving radius (-Rasy) of a blastwave solution in a flat universe dominated 
by collisionless matter depends upon Eq, Qigm, zi, Ho, and also on the form of Rs{t). When 
the other parameters are equal, solutions that expand more slowly reach larger values of Rgsy 
This is why the V95 and OM88 solutions eventually differ. To sec how the shock expansion rate 
affects i?asy) consider how the initial energy Eq shifts from the IGM to the collisionless component, 
assuming that the shock sweeps the IGM into a thin shell. As an IGM shell of mass M overtakes 
a collisionless shell of mass dM^, the gravitational energy of the collisionless shell increases by 
GM dMd/Rg and the gravitational energy of the IGM shell decreases by the same amount. In 
terms of r){x) = Rg/Rasy and x = {1 + z)/{l + zi), we can write the energy transfer rate as 

dE _ 3H^Rl^n,GM{l + z,) 

dx~ 4G ■ ^ ' 

The total amount of energy transferred approaches Eq &s x ^ 0; hence, 

AGEq 



-f^asy — 







.3i/o^OlGM(l + ^l 

From § 2.1 we have r]{x) = (1 — x^/^)^/^ when VLq = 1, which leads to 



1/5 / /.I X -1/5 

xTj^T]' dx . (2-23) 



AOGEo 



H^niGuil + zi) 



1/5 



(2-24) 



This is identical to the result obtained by Voit (1994) in the thin-shell approximation. The value 
of Rgsy from § 2.1 is 11% larger than this because the IGM density in the exact solution falls 
off more gradually behind the shock front, slightly reducing the amount of gravitational energy 



associated with a given asymptotic radius. Note that blastwaves starting at later times expand to 
larger comoving radii. Ostriker k, McKee (1988) took 



■n(x)= (l-0.66x2/2)(l-x^/2) 



2] 1/5 



(2-25) 



which yields an asymptotic radius in the thin-shell approximation 16% smaller than that in 
equation (27). The actual asymptotic radius of the OM88 solution is 20% smaller than the value 
of i?asy from § 2.1. Figure 3 shows that the OM88 form for r/(x) increases more rapidly at early 
times, transferring the initial energy to the dark matter sooner, when the energy needed to achieve 
a given comoving radius is larger. 



Although the blastwave radii of the approximate solutions do not differ dramatically, the 
predicted shock velocities can be quite different at late times. Figure 4 displays the shock velocities 
corresponding to the approximations illustrated in Figure 3. Because the OM88 approximation 
approaches its asymptotic radius so quickly, the shock velocity drops off rapidly beyond t ~ Ati. 
In this approximation the shock cools prematurely. The VB87 approximation is more accurate 
and remains close to V95 until t ~ 30ti (~ ^ij^j^^ti for Q-igm = 0.1.) 



Immediately behind the shock front the gas density is four times the current IGM density. 
Because the solution of § 2.1 is directly analogous to the Scdov solution, the postshock density 
drops as r decreases. Gravitational forces can compress the shell further, but the compression 
does not become significant until t ~ f^jQM ii (§ 2.3). Ikeuchi et al. (1981) argue that explosions 
in the IGM should eventually produce dense shells, but postshock densities in their numerical 
calculations do not exceed 4/9igm until t ^ lOti . 



2.5.2. Shock Velocity 



2.5.3. Postshock Density 



3. 



Ionization 



At z < 5, the intergalactic medium appears to be highly ionized (Gunn k, Peterson 1965; 
Steidel & Sargent 1987; Webb et al. 1992; Giallongo et al. 1994). Blastwaves at these redshifts 
can expand adiabatically until radiative cooling limits their growth (§ 2.2). When the IGM is not 
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yet fully ionized, shocks initially propagate into a cosmological H II region created by ionizing 

radiation from the source of the blastwave, prcsTimably a quasar or a protogalactic starburst. If the 
blastwave energy exceeds the ionizing output of the source, the shock can potentiahy enlarge the 
ionized cavity. Shocks reaching the neutral IGM at speeds greater than ~ 50kms~^ coUisionally 
ionize additional gas. This section determines the conditions under which shock heating adds to 
the total mass of ionized gas. 

Ultraviolet radiation from a quasar or starburst in the early universe quickly ionizes the 
surrounding IGM. Analytical models for the expansion of cosmological H II regions have been 
developed by Donahue & Shull (1987), Shapiro & Giroux (1987), and Madau & Meiksin (1991). 
The recombination time in a uniform IGM is (6.3 x 10i°yr)(l + zy^fl^^ly^h^^T^-^, where 
T4 = T/10^ K. Thus, at z < 7, ionized intergalactic gas with I^igm < 0.1 remains predominantly 
ionized, and a burst of ionizing radiation produces a permanent H II region. Although the neutral 
fraction in this region does not necessarily remain small enough to satisfy existing constraints on 
the total optical depth of intergalactic Lya, it does permit blastwaves to propagate at speeds slower 
than 80kms~^ without suffering heavy ionization losses. If the central source emits (10^^ erg)^76i 
in ionizing photons of average energy (13.6eV)£'R,yd, the total gas mass of the ionized region is 



Mi = (5.4 X 10^^ Mq) Eje^E^^^ , 



(3-1) 



assuming primordial abundances. This mass corresponds to a comoving radius 
^ {12.3Mpc)e]^JIe^^^^ H^q^^Q^qI^ . The H II region grows much more rapidly than any 
hydrodynamic disturbance until the ionizing source shuts off and grows slowly thereafter as diffuse 
recombination radiation continues to provide ionizing photons. After the ionizing episode, the 
warm ionized gas can drive a 10 — 20kms~^ shock into the neighboring neutral IGM. Madau & 
Meiksin (1991) employed a special case of the transformation described in § 2.1 in their analytical 
models of these shocks. 

CoUisional ionization adds to the ionized mass when the blastwave energy exceeds the 
radiative ionizing energy and the energy of the blastwave is relatively small. The mass of shocked 
gas inside an adiabatic cosmological blastwave is Masy?7^, where 



2G 



= (1.2x10^^ Me 



e) 



1 + Oo-^l 
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in a flat universe. Thus, an adiabatic blastwave eventually expands beyond the preionized region if 



Eqi 



< 0.22£jRy(J 



Eeih^oQ 



IGM 



(1 + noZi)3/2 
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(3-3) 



Blastwaves generated by large-scale energy sources {Eqi ^ 1) remain within the preionized region 
unless the initial burst of ionizing radiation is quite weak. Smaller-scale blastwaves can escape 
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the preionized region more easily. The thermal energy density immediately behind a ~ 50kms~ 
shock is just enough to ionize all the incoming hydrogen. Below this threshold, the postshock gas 
remains mostly neutral. Let ry/ be the dimensionless comoving radius at which Vg = 50kms~^, so 
that Mi^y-qj is the total mass of gas a cosmological blastwave can ionize. Figures 5 and 6 show how 
7]! and Masyr?f depend on the parameter combination £"61/150(1 + -2i)^^^^igm ^ universe. 
Blastwaves in an open universe behave similarly. 

When 77/ <C 1 , a large proportion of the blastwave's energy can go into ionization of hydrogen. 
Small-scale blastwaves propagating into a neutral IGM lose significant energy to ionization when 
t^s ~ 50 — SOkms"-*^ and cool via Lya when t>s ~ 10 — SOkms"-*^, as long as JIigm <; 0.01. The 
blastwave switches from an adiabatic solution to a pressure-driven snowplow solution at r] = rji, 
so the total mass within the blastwave asymptotically approaches Masy??/^^ (§ 2.2). In blastwaves 
with rji approaching unity, most of the initial energy goes into the collisionless component when 
Oo ~ 1 or into kinetic energy that redshifts away when ^ 1, leaving little available for 
ionization. Consequently, large-scale explosions are quite inefficient at ionizing the IGM. 

Around quasars, photoionization exceeds collisional ionization. Although some quasars might 
expel winds with luminosities comparable to their radiative luminosities (Voit, Weymann, &; 
Korista 1993), constraints on the distortion of the microwave background spectrum imply that the 
total energy in quasar blastwaves cannot exceed the total radiative output of quasars by a large 
factor (Voit 1994). Since the integrated radiative output of a typical quasar is ~ 10^"^ erg, quasar 
blastwaves remain within their preionized cavities (see Fig. 6). 

The ionizing properties of protogalaxies are harder to constrain because it is difficult to 
measure the fraction of ionizing photons that escape. In a starburst, the ratio of hydrodynamic 
to photoionizing output, integrated over the burst, can be as large as 0.3 (Leitherer & Heckman 
1995). If the galaxy itself absorbs > 70% of the ionizing photons while allowing all of the 
hydrodynamic energy to escape, collisional ionization could, in principle, exceed photoionization 
in the surrounding IGM. For this to happen, the total magnitude of the protogalactic starburst 
would have to be quite small. If r^icM =0.1, hr^Q = 1, and Ejqi/ E^y^^EQi = 0.3, a blastwave 
beginning at zi = 5 would not augment the photoionized cavity unless Eqi < 10~*^'^, implying 
< 10^ supernovae. 



4. Comments on the Fluid Approximation 

Most investigations of intergalactic hydrodynamics, including this one, assume that the fluid 
equations adequately describe the behavior of the IGM, yet the mean free paths of fast particles 
in the IGM are not necessarily small. Coulomb collisions with electrons in a 10^ K plasma stop 
a proton moving at lO^ugkms"^ over an electron column density ~ lO^^-ygCm"^ (Spitzer 1962). 
This column corresponds to a comoving distance ~ (0.1 Mpc)(J7igm/10~^)~^V8. In a hot plasma 
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at 10 T7 K, the column a proton traverses before Coulomb collisions significantly deflect it is 
~ 10^^ T7 cm^^ (Spitzer 1962), corresponding to a comoving distance ~ (1 Mpc)(r2iGM/10~^)~^T7 . 
These distances are similar to the length scales that emerge from cosmological blastwave solutions. 

An intergalactic magnetic field of lO^^*' B20 G confines protons moving at 10^ vg K to within 
a gyroradius of 0.3 fs-B^^ Mpc, so a primordial field exceeding 10~^^ G would restore the fluid 
approximation on intergalactic scales. The actual strength of the intergalactic magnetic field is 
hard to constrain (Kronberg 1994). Because the current Galactic field (~ 10~^ G) must arise 
from an initially small seed field, some investigators have evaluated the primordial field strengths 
required by various models for Galactic field generation. Microwave background scattering inhibits 
the motions of electrons in rotating protogalaxies, causing the electrons to lag behind the ions. 
This weak current generates a small (10~^^ G) intergalactic field (Mishustin &: Ruzmaikin 1972). 
Galactic dynamo action can potentially amplify a seed field this small by a factor ~ e^''°' ~ 10^^, 
where A^rot ~ 30 is the number of times the Galaxy has rotated (Parker 1979). Kulsrud & 
Anderson (1992) have criticized this picture, arguing that is it very difficult for a galactic dynamo 
to create the current field from a very small seed field (but see Field 1995). The Galactic field 
would then have to arise from a much larger (~ 10~^^ G) mean primordial field, amplifed primarily 
by the initial compression of the protogalactic gasj^ Uncertain detections of Faraday rotation in 
radio-selected QSOs behind damped Lya systems at z ^ 2 suggest that early galactic disks had 
fields ~ 10~® G when Afot <C 30 and would support the view that the intergalactic field must be 
» 10-21 G (Wolfe, Lanzetta, & Oren 1992). 

In interstellar and interplanetary shocks, magnetic fields and plasma turbulence transfer 
energy between particles more efficiently than Coulomb collisions do, allowing postshock fluid 
structures to develop on length scales much shorter than the Coulomb mean free path (Draine &: 
McKee 1993). Even in the absence of a magnetic field, electrostatic shocks mediated by plasma 
turbulence can thermalize the incoming particle energy in a distance similar to the Mach number 
times the Debye length (Tidman &; Krall 1971). Since such shocks are quite thin, even in the 
IGM, the radii of cosmological blastwaves are still well defined and we can safely use the fluid 
equations and strong-shock jump conditions to describe how they propagate. 



5. Global Effects 



At redshift z ~ 4 the intergalactic medium is 
clouds have formed. If cosmological blastwaves are 
of the IGM or for producing Lya clouds, they must 
electrons inside such blastwaves scatter microwave 



already highly ionized, and Lya- absorbing 
responsible for maintaining the high ionization 
have filled intergalactic space by this time. Hot 
background photons, distorting the microwave 



Certain scenarios for inflation in the early universe can produce primordial fields as high as 10 ^ G (Ratra 1992). 
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background spectrum (Weymann 1965, 1966; Sunyaev & Zeldovich 1980). Limits on this distortion 
restrict the total blastwave energy density in the early universe. Metallicity limits further constrain 
the energy density in supernova blastwaves, if the metal abundances of Lya clouds reflect the 
average metallicity of the IGM. The metallicity constraint allows adiabatic blastwaves beginning 
at zi < 7 to shock most of the IGM by z ~ 4 only if the energy per blastwave is small and the 
comoving number density of the blastwave sources is high (~ IMpc^^). Thus, the metals in Lya 
clouds were probably produced in small-scale starbursts (< 10^^ erg) rather than in large-scale 
protogalactic superwinds (> 10^'' erg). This section outlines these constraints in more detail. 



5.1. Comptonization 

The COBE satellite has shown that the cosmic microwave background spectrum deviates 
little from a pure blackbody. This lack of deviation can be used to limit the amount of energy 
injected into the IGM by cosmological blastwaves. Stars with a standard initial mass function 
release a fraction esN ~ 5 x 10~^ of their initial rest mass energy in the form of kinetic energy 
from supernovae. The mean thermal energy density introduced into the IGM by supernovae in a 
star forming episode at zi is then 

Uiizi) = (1.5 X 10~^^ergcm-^){l + zifhlj,A , (5-1) 

where /esc is the fraction of the supernova energy that escapes a star-forming galaxy and $7* is 
the mass density of stars formed at zi, in units of the critical density.0 Since Eq is constant in 
adiabatic cosmological blastwaves, Ui{z) = ?7i(zi)[(l + z)/{l + ^i)]^. Shocks in a low-density, 
preionized IGM remain adiabatic until z has changed significantly (Figure 2). The amount of 
y-distortion expected from such a star-forming episode is 

yi « (2.2 X 10-5) hofesA[{l + zif/^ - (1 + zi)-2] , (5-2) 

assuming n = 1. COBE finds y < 2.5 x 10"^ at the 2a level (Mather et al. 1994), implying 

2/5 

^*h5ofesc < 0.05 - 0.08 for zi = 7 - 5. At zi > Shl'^ - 1 in a flat universe, Compton 
scattering transfers all the supernova energy to the microwave background, resulting in the limit 
17,/csc< 2.8 X 10-3 (1 + zi). 

Fluctuations in the numbers of blastwaves along different lines of sight will cause the 
brightness of the microwave background to vary on small angular scales (Hogan 1984; Yoshioka 
& Ikeuchi 1987, 1988), but these fluctuations are relatively insignificant as long as the blastwave 
energies do not greatly exceed 10^^ erg. A single blastwave distorts the microwave background by 



^ About 72% of the initial energy passes into thermal energy as the blastwave develops (Sedov 1959). 
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which reduces the background's brightness temperature by AT/T ~ —2Ay. The temperature 
fluctuations induced by individual blastwaves are unobservable unless r] is very small or Eqi is 
very large. 

Brightness temperature fluctuations from an ensemble of blastwaves depend upon how many 
blasts lie along a typical line of sight. The cross section of a blastwave is oc ?7^, but the amount 
of time a blastwave spends at a size 77 ^ 1 is oc 77^/^. The number of blastwaves of size 77 
along a line of sight is therefore J\fjj oc rj^/^. Large blastwaves contribute most heavily to y, so 
A/'r,(r7 ~ 1) ^ y/Ay and 



ATr, < (3 X 10^ 



»g(l + ^i)' 

1 -I- fioZi 
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When there are many blasts along a line of sight, the y fluctuations from blasts of size 7/ are 
~ A/"j|/^Ay oc 77^/^. Blasts with 77 ~ 1 thus dominate the temperature fluctuations, which are of 
order 



AT 



< (3 X 10" 



-7\ 



nl/5 



l + Z 
l + Zl 



p3/10, 4/5^1/5 
-^61 "50 ^'IGM 



1 -I- ^qZi 

The temperature fluctuations induced by galaxy-scale explosions, which cannot exceed 
jQ]^ Lui ^1 — 5 and JIigm ~ 0.01, are presently unobservable. 



(5-5) 



10-7^;^/^° for zi 



Brightness temperature fluctuations imprinted by individual young (77 <C 1) blastwaves can 
exceed AT/T ~ 10"^ but these are rare. When f^iGM ~ 0.01, explosions with 77 < 0.02£'g(^° can 
produce temperature fluctuations at this level. The angular sizes of such blasts at 2 ~ 4 in a flat 



1 /2 ' ^ 

universe would be ~ 5 Eq[ arcsec, and they would cover a fraction < (3 x 10-'^)£;6i of the sky. 



.3/4 



5.2. Metallicity 



Supernova-driven winds from galaxies carry metals into intergalactic space, and the fraction 
of the metal-bearing supernova ejecta that escapes a starbursting galaxy should be ~ /esc- Since 
stars formed with a normal initial mass function eventually expel about 1% of the input mass in 
the form of metals, a star- formation episode enriches the IGM to a mean level ~ /esc^*/^iGM 
times solar metallicity. If the enrichments of Lya clouds reflect average IGM abundances, then 
upper limits on the intergalactic metallicity constrain the energy input from supernova-driven 
blastwaves. Lu (1991) derived Lya-cloud metallicitics ~ 10-^ times solar from stacked C IV 
spectra of clouds with neutral columns exceeding A'hi > 10^^ cm-^. Recent Keck spectra of 
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smaller clouds (Cowie et al. 1995) now indicate metallicities 10 ^ — 10 times solar. Similar 
metallicities in the global IGM would imply /esc^^*/^iGM ^ 10"^. 

This constraint assumes that metals escaping from protogalaxies are evenly mixed into the 
IGM and Lya clouds. The assumption of even mixing is most justifiable if protogalactic blastwaves 
filled the IGM before Lya clouds condensed. If mctal-bcaring blasts encompassed only a small 
fraction of the IGM by the time Lya clouds formed, the variance in intergalactic metallicities 
would be high. Enriched clouds would have high metallicities, while the others would be metal-free. 
Tytler &; Fan (1994) determined that Lya clouds with metallicities 10~^ times solar must be 
extremely rare, and Cowie et al. (1995) find a surprising consistency in cloud metallicities. While 
it is possible that the clouds are self-enriched by associated star formation (e.g. Cowie et al. 1995; 
Madau & Shull 1995) and do not reflect the metallicity of the IGM, they would still provide an 
upper bound on IGM enrichment at the time of cloud formation. The metallicity constraint would 
not apply if Lya clouds formed before protogalactic blastwaves reheated the IGM and avoided 
mixing with the hotter metal-rich gas as it swept past. Nevertheless, further explorations of 
metallicity variations in Lya clouds might provide very useful constraints on the star-formation 
history of the universe. 



5.3. Porosity 



Blastwaves initiated after z = 7 can fill the IGM by z = 4, while still satisfying the above 
constraints on microwave background distortion and metallicity, as long as the typical explosion 
size is relatively small. Blastwaves of characteristic energy scale Eq stemming from star formation 
at zi have a comoving number density = (6.3 x 10"^ Mpc^^)r2*/esc-^!'6i- Multiplying this number 
density by the typical blastwave volume gives a porosity parameter 

Q = 111 n,UcEei^'h^o^'n-^{^{l + zi)-3/5^3 (5_Q^ 

analogous to the porosity parameter for supernova remnants in the Galactic interstellar medium 
(Cox & Smith 1974). Blastwave sources distributed randomly in space shock a fraction ~ 1 — e^'^ 
of the total IGM volume. Shocks from clustered sources merge early, increasing the effective Eq 
and lowering the porosity, so Q can be considered an upper limit on the true porosity. Writing Q 
in terms of the typical blastwave shock velocity Vs yields a particularly illuminating expression: 

where x = (1 + z)/(l + zi). The initial energy input decays adiabatically (oc x^), and Q approaches 
unity when this energy has been distributed evenly throughout the IGM. The Comptonization 

^ 1/2 3/2 

limit on /esc^* at zi = 7 permits Q > 1 when Vg < (180kms~ )r2jQj^ h^Q x. The more stringent 
metallicity limit permits Q > 1 when Vg < (79kms~-'^)x. 
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More generally, let Z^i be the IGM metallicity in units of 10~^ times the solar metallicity. 

Spreading the blastwave energy throughout the IGM gives an average IGM temperature 
~ (2 X 10^K)Z_2a;^, if there arc no ionization or radiative losses. Requiring Q > 1 at 
z = 4 for z\ = 7 in a flat universe then implies E'o < (4 x 10^'^ erg)(5~^/^Z^2^/igQ^r2iGM and 
rift > (1.6 Mpc^^)/i50'2^-2'^^- This density is much larger than the number density of galaxies 
and is comparable to the number density of Lya clouds (§ 6.1). Figures 7 and 8 illustrate these 
constraints graphically for fiioM = 0.01 and 0.03. The enrichment levels of Lya clouds thus 
suggest that the metals they contain originated in relatively small-scale starbursts close to the 
clouds themselves or in individual population 3 supernovae rather than in enormous protogalactic 
starbursts. 

Quasar blastwavcs do not necessarily enrich the IGM with metals, thus they are not bound 
by the metallicity limit. The Comptonization limit, however, still applies. Because the number 
density of quasars is low, the energy scale of their blastwaves must be high, if these blastwaves 
are to fill the universe quickly. The energy density needed for quasar blastwaves to give Q > 1 by 
z ~ 4 violates the COBE limits on y-distortion (Voit 1994; see also Figures 7 and 8). Apparently, 
quasars are too rare to shock-heat the entire IGM before Lya clouds appear. 



6. Blastwaves and Lya Clouds 



Although we have been assuming the IGM is homogeneous, intergalactic space is obviously 
inhomogeneous. The IGM is filled with clouds responsible for the numerous Lya absorption 
lines, known collectively as the Lya forest, that litter the spectra of high-redshift quasars. The 
Lya clouds remain puzzling. Many authors have speculated that these clouds might be linked 
somehow to cosmological blastwaves. Some have suggested that the absorption lines come from 
the intergalactic shock fronts themselves (e.g. Chernomordik &; Ozernoy 1983). They could 
also be condensed fragments of primordial blastwaves, now pressure-confined by the IGM (e.g. 
Ostriker h Cowie 1981; Ikeuchi h Ostriker 1986; Vishniac & Bust 1987), gas that has settled 
into shallow mini-halos of dark matter (Rees 1986; Bond, Szalay, & Silk 1988), gas collecting at 
caustics in the dark-matter velocity field (McGill 1990), or pancakes of shocked gas where that 
velocity field converges (Cen et al. 1994). This section explores whether cosmological blastwaves 
produce observable Lya absorption lines and examines how these shocks interact with preexisting 
intergalactic clouds. We find that the postshock H I column densities (A^^hi) in the ambient 
IGM are much smaller than those of intergalactic Lya clouds, but blastwaves that interact with 
preexisting clouds can amplify the Lya equivalent widths of these clouds significantly. 



6.1. Lya Cloud Characteristics 
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The joint distribution function describing the number of Lya clouds per redshift interval 
per column density interval can be written 

M^N^^H + zV , (6-1) 



dN-aidz 

where = Nm/W^^ cm.~'^ . Marginally opaque clouds ~ 1) absorb most of the photons 
when 1 < /5 < 2. For this reason, clouds with A'^hi ~ 10^^ cm~^ are considered typical, even though 
a power law distribution of column densities fits the data. Press & Rybicki (1993) find TVq « 2.6 
and P = 1.43 it 0.04 in a column density range at least as large as 10^^ cm~^ ^ A^hi ^ 10^^ cm~^ 
for 7 = 2.46 at redshifts 2.5 ^ z ^ 4 (Press, Rybicki, & Schneider 1993). The average comovmg 
distance between clouds with A^hi > 10^"^ cm^^ is then ~ (7Mpc)[(l + 2.5fio)/3.5]"i/2/^5o^ at 
z = 2.5 and ~ (1.7Mpc)[(l + 40o)/5]"^/^/i5o^ at z = 4. 

Near the redshift of the background quasar, the number density of Lya absorption lines 
decreases (Weymann, Carswell, & Smith 1981; Murdoch et al. 1986). This "proximity effect" 
presumably occurs where the ionizing radiation field of the quasar overwhelms the ambient 
ionizing field. The mean intensity of ionizing radiation at 2: ~ 3 inferred from the proximity effect 
is (10~^^ ergcm"^ s~^ Hz~^ ster~^) J21, where J21 ~ 1 — 3 (Bajtlik, Duncan, & Ostriker 1988; 
Bechtold 1994). In photoionization equilibrium, the fractional abundance of atomic hydrogen in 
an intergalactic cloud is 

/ho « 2.8 X 10-^ (1 + zf D,iniGMhloJ2iV-^ , (6-2) 

where -Del is the ratio of the cloud density to the IGM density and T4 is the cloud temperature in 
units of 10^ K, giving a total column density for a typical cloud of 

A^H = (3.6 X 10^0 cm-2) (1 + z)-^ NuD-.^^^^^h^^T^'^ J21 . (6-3) 

Recent observations of the quasar pair Q1343-I-2640A, B, indicate that the transverse cloud radius 
Tci = (200 kpc)/i^Q^r2oo for clouds with A''i4 ~ 1 lies in the range 80/150^ kpc < < 560 /i^g^ kpc at 
z ~ 1.8 (Bechtold et al. 1994; see also Dinshaw et al. 1994). This result implies 

D,, ~ 0.8<'^GV/i-o'/'Ti'-V2YV2"of (6-4) 
ATh ~ (2 X 10^9 cm-2)</^/.-o^/^rr 4/^2^0' (6-5) 

at 2; ~ 2, if the clouds are roughly spherical. Sheetlike clouds would have a larger D^i and 
smaller A'h- Assuming spherical clouds, Bechtold et al. (1994) derive a comoving number density 
~ (2 Mpc~'^)/i5o'"20C) fo^ clouds of A'^14 ~ 1 at z ~ 2. Perhaps coincidentally, this density is similar 
to the comoving number density of blastwave sources needed to enrich the entire IGM to 10~^ 
times solar by this redshift (§ 5.3). 



^The redshift evolution parameter 7 is still not well constrained; current values range from 1.89 (Bechtold 1994) 
to 2.75 (Lu, Wolfe, & Turnshek 1991). 
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6.2. IGM Ionization Constraints 



We see a Lya forest because the spaces between Lya clouds are either gas free or highly 
ionized. Gunn 8z Peterson (1965) inferred /o^^igm = ^^H i < 6 x 10~^ h'^^ from the lack of 
continuous absorption at z ~ 2. Observations at the wavelength of He II Lya (304 A) toward 
a quasar at z = 3.286 show that intergalactic He II absorbing gas fills velocity space a,t z ^ A 
(Jakobsen et al. 1994). Although this observation appears to confirm that the IGM has a spatially 
continuous component, the cumulative He II absorption of many very thin velocity-broadened 
Lya clouds (TVhi ;^ 10~^^ cm~^) could also produce the observed signal (Madau k, Meiksin 1994). 
Current limits on JIhi = 10~^ fig iii ^ spatially continuous component arc Vl^ < A.Sh'^Q at z ~ 2.5 
(Steidel & Sargent 1987) and < 2.5 h^Q at z ^ A (Webb et al. 1992; Giallongo et al. 1994). 
Photoionization equilibrium thus implies JIigm < 0.016 /i5Q^r4-^J2/^rig^^ at z = 4.1. 

Alternatively, the IGM could be collisionally ionized by early supernovae (Tegmark, Silk, & 
Evrard 1993), but if the metallicity of the IGM is similar to that of Lya clouds, photoionization 
is likely to be more important. The mean temperature in an IGM heated by supernova-driven 
blastwaves is ^ (2 x 10^K)Z_2 (§ 5.3), and pure collisional ionization at 2 x 10^ K requires 
^^IGM < 3.6 X 10~^ ris for Lya to be transparent at the observed levels, a condition more stringent 
than the photoionization condition. An elevated temperature does, however, reduce the IGM 
recombination rate, easing the photoionization constraint on flicM (e.g. Shapiro 1995). 
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6.3. Shocked Shells as Lya Clouds 

In principle, a sufficiently large compression of the IGM in a shock along the line of sight to a 
quasar could produce an observable Lya feature (Ozernoy & Chernomordik 1978; Chernomordik 
& Ozernoy 1983). This can happen if f^iGM ~ !> but since f^iGM ^ 10^^, as the ionization 
constraints indicate, the H 1 column densities of cosmological blastwaves in the ambient IGM are 
much smaller than typical Lya-cloud column densities (A^^hi ~ lO^^cm"^). The total hydrogen 
column density within a cosmological shock is Nasy{l + z)'^r), where 

TVasy = (1.6 X W'^ cm-')Eli''h%''nl[;i,{l + z,f/' . (6-6) 

Figure 9 shows how the shocked column density evolves for various blastwave parameters 

—3/2 

and r^iGM = 0.01. Behind a strong adiabatic cosmological shock, D^i = 4 until t ~ ^igm 
(§ 2.5.3). Initially, electron collisions keep the postshock gas highly ionized. As the shock 
slows, photoionization maintains a relatively high ionization level. If collisional ionization can be 
neglected, the neutral column behind the shock is 

AThi « (1.7 X 10^^ cm-2)(l + zfEli'hlTn',l^^T,-'-'j^,\l + z^fl^ ■ (6-7) 

Figure 10 illustrates A^hi(-z) for a range of blastwave parameters and Oigm = 0.01. The lines 
indicating AThi are dashed when collisional ionization exceeds photoionization and solid when 
photoionization exceeds collisional ionization. Low-energy blastwaves produce larger neutral 
columns because they have smaller postshock temperatures (T4 oc -Eg(^). Unless J^igm ~ 1; Lya 
absorption in adiabatic blastwaves is much weaker than in Lya clouds. 

Cooling can potentially magnify the postshock density, boosting the neutral column in the 
blastwave, but the compression ratios of radiative cosmological shocks in the ambient IGM are 
unlikely to be very large. At z 2, shocks with E^) <^ 10^"^ erg can cool in a Hubble time in a 
medium with Ojgm 0.03 (Figure 2). The requisite shock velocities are ;^ 40 — 50kms^^, not 
much larger than the velocity dispersions of Lya clouds. The widths of intergalactic Lya lines, 
if thermal, imply temperatures ^ 3 x 10^ K, consistent with photoelectric heating by quasars. If 
the intercloud gas is at least as hot, radiative IGM shocks have Mach numbers <^ 3, too weak to 
produce iVni ~ 10^^ cm~^ when r^iGM ^ 0.03. Allowing I^igm ^ 0.1 requires an IGM temperature 
^ 10^ K for Lya to remain transparent. Even in this denser environment, the Mach numbers 
of radiative shocks, with velocities < 200 km s^^ (§ 2.2), would still be too low to increase 
sufficiently. 



6.4. Shock-Cloud Interactions 

Since blastwaves in a low-density IGM have insignificant neutral columns, Lya clouds form 
through some other means. The porosity and metallicity constraints discussed in § 5 argue 
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against pressure confinement of Lya clouds. Mechanisms in which interstellar gas collects in 

"minihaloes" of dark matter (e.g. Rees 1986) or converges at caustics in the large-scale velocity 
field (e.g. McGill 1990, Cen et al. 1994) appear to be more likely. Fluctuations in the spatial 
distribution of collisionless matter certainly alter the density field of the IGM, which we have so 
far idealized as uniform. Real intergalactic shocks propagate through an inhomogeneous medium; 
hence, the solutions derived above become accurate only after the shock radius exceeds the length 
scales of non-linear density perturbations. The remainder of this section briefly examines how 
interactions between cosmological blastwaves and density inhomogeneities proceed and suggests 
that intergalactic blastwaves might amplify the neutral column densities of preexisting Lya clouds. 

Shock propagation into an inhomogeneous interstellar medium has been analyzed extensively 
(McKee k Cowie 1975; Woodward 1976; Nittman, Falle, & Gaskcfi 1982; Stone & Norman 1992; 
Klein, McKee, & Colella 1994). When a shock of velocity Vs encounters a cloud -Del times denser 

— 1/2 

than the ambient medium, it drives a shock of velocity D^j Vs into the cloud. In a cloud-crushing 
time {tec ~ -Cy^rciWci), the internal shock crosses the cloud. Detailed numerical simulations show 
that the vorticity generated in the cloud-shock collision shreds the cloud in a few cloud-crushing 
times (Nittman et al. 1982; Stone & Norman 1992; Klein et al. 1994). If Lya clouds are indeed 
several hundred kpc in size, the cloud shredding time approaches a Hubble time for Vg less than a 
few times lO^kms"^. Since the postshock cooling time in the cloud is smaller than the ambient 
cooling time by a factor ~ D^j, such long-lived clouds, or fragments thereof, might collapse to 
much higher density before vorticity destroys the cloud. The resulting AThi would then increase 
by a factor ~ Vg/v^f^, where Vth ~ 30kms~^ is the thermal velocity dispersion of the photoionized 
cloud after it has cooled to its equilibrium temperature. Such a shock would impart a velocity 

— 1/2 

D^^ Vg to the cloud, perhaps carrying it out of its potential well. Figure 11 illustrates how large 
-Del must be for shocked clouds to cool within a Hubble time. In the example shown, shocks of 
Eq = 10^^, 10^^, and 10^^ erg begin propagating at = 4 in an Q = 1, /150 = 1 universe with 
^^IGM = 0.01. Clouds with -D^ = 6 at 2; ^ 2, for instance, can cool after encountering shocks of 
Vs ^ 200 km s-^ 

Computing in detail how cosmological blastwaves interact with developing gaseous structures 
in the IGM is beyond the scope of this paper; however, the cooling times estimated here indicate 
that blastwaves could dramatically alter the observable characteristics of preexisting Lya clouds. 
In essence, shocked clouds are pressure-confined, but the pressurizing medium does not pervade 
the universe. Instead, the pressurized regions are the interiors of individual blastwaves. Clouds 
near blastwave sources would then have systematically higher neutral column densities than clouds 
in the ambient IGM. 



7. Summary 
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This paper has investigated the behavior of protogalactic and quasar-driven blastwaves in an 
expanding universe with J^igm Oq. It presents a new analytical solution for adiabatic blastwaves 

—3/2 

that applies until t ~ tiiljQj^j , when self-gravity starts to compress the postshock gas into a thin 
shell. Although two different timescales {t and ti) enter into the problem, a self-similar solution 
still exists. The similarity of the universe's expansion allows us to define another time coordinate 
(t) in which the Sedov- Taylor solution satisfies an appropriately transformed set of fluid equations 
if 7 = 5/3. Cosmological blastwaves thus begin as Sedov- Taylor blastwaves and approach a 
constant comoving radius as t grows large. While t <C tiQ^^^^ , the gas density immediately behind 
the shock remains equal to four times the ambient density, so cosmological shocks do not cool 
until they are moving quite slowly. 

Applying this solution to protogalactic starbursts and quasars shows that small-scale 
explosions 10^'' erg) are probably more important than large-scale explosions 10^^ erg) in 
determining the state of the IGM and its Lya clouds. Large-scale explosions are quite inefficient at 
ionizing the IGM. Cosmological effects rob most of the energy from large-scale blastwaves before 
these shocks slow to speeds that ionize efficiently lOOkms"^). In addition, since the comoving 

—2/5 

volume per unit energy of a blastwave is cx ; small-scale explosions can fill the IGM much 

more easily. Quasar-driven blastwaves are too rare to fill the IGM by z 4 without violating the 
y-parameter limit from COBE. Large-scale protogalactic blastwaves that filled the universe by 
such an early time would raise the mean metallicity of the IGM well above the values observed in 
Lya clouds at z ~ 2 — 3. 

Contrary to some suggestions that Lya clouds might be shells of gas shocked by cosmological 
blastwaves, the solution derived here shows that the postshock gas never becomes dense enough 
to produce iVni <; 10^^ cm"^ when Oigm ^ 0.03. Even if the blastwave has slowed enough to cool, 

the compression of the ambient IGM is still too small to raise A^hi to the levels observed in Lya 
clouds. Blasts that pass through over dense regions of the IGM, such as preexisting Lya clouds, 
might still produce observable signatures. Compression and cooling can dramatically amplify the 
neutral fraction in a shocked cloud, so clouds engulfed by blastwaves could have systematically 
higher values of Nm- This effect would increase the probability of finding clouds near blastwave 
sources. 

During the development of this paper, conversations with Steve Balbus, Megan Donahue, 
Mike Fall, Piero Madau, and Sterl Phinney were particularly helpful. The author gratefully 
acknowledges the hospitality of the Virginia Institute for Theoretical Astrophysics during the early 
stages of this project. The work presented here was supported by NASA through grant number 
HF-1054. 01-93 from the Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 1. — Growth of an adiabatic cosmological blastwave when JIigm ^ f^o = 1- This figure 

illustrates how quickly the radius and mass of an adiabatic cosmological blastwave arising from 
an intergalactic explosion at cosmic time ti approach their asymptotic values. At cosmic time t, 
the shock radius has reached a fraction i?/i?asy of its asymptotic comoving value and contains a 
fraction M/M^y of the mass it will ultimately encompass. 
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Fig. 2. — Shock velocities of adiabatic cosmological blastwaves. The velocity Vs(z) of the adiabatic 
shock driven into the IGM when an energy (10^^ eigjE^i is introduced at redshift zi depends on zi, 
O, and the value of E'ei/iso^^iGM- Here we show how Vg varies with z for several values of £^61^^igm 
when /i5o = 1,0 = 1, and r^ioM ^ 1- Curves corresponding to shocks beginning at three different 
initial redshifts (^i = 1, 4, 7) are shown. The solid lines indicate where cooling begins to violate the 
assumption that the shocks are adiabatic. The vertical line at z = 7 shows where Compton cooling 
against the microwave background becomes important, and the two lines across the bottom show 
where radiative cooling becomes important for Ojgm = 0.01 and 0.03, assuming an IGM metallicity 
of 0.01 times solar (radiative cooling functions from Sutherland &; Dopita 1993). The horizontal 
dashed line shows where ionization losses would become important if the ambient medium were 
neutral. Since the Vs{z) curves lie primarily above and to the right of the solid lines, cosmological 
shocks at z < 7 are adiabatic for a wide range of parameters. 
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Fig. 3. — Approximate solutions for adiabatic cosmological blastwaves when JIigm <IC fio = 1- Solid 
lines (V95) show the exact solutions from § 2.1 for Rs{t) in comoving coordinates, given a 10®^ erg 
blastwave starting at = 5 in intergalactic media with Oigm = 0.1 and fiicM = 0.01. These 
solutions become invalid when t ^ h^iQu- Dotted lines (OM88) show the approximate solutions 
derived by Ostriker & McKee (1988) for these same parameters. The long-dashed line (VB87) gives 
the approximate fit of Vishniac & Bust (1987) to numerical models. The short-dashed line (B85) 

3/2 

gives the self-similar solution from Bertschinger (1985), valid when t ^ ^i^igm shocked 
IGM comoves with a coUisionless shell that has responded to the blastwave. 
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Fig. 4. — Approximate shock velocities. The hnes in this figure trace the shock velocities given by 
the approximate blastwave solutions in Figure 3. 
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Fig. 5. — Dimensionless comoving radius at which shock ionization becomes ineffective {Q = 1). 
The quantity 77/ is the value of 77 = Rg/Rasy at which the shock velocity drops below 50kms~^. 
Beyond this point, the shock no longer provides enough energy to ionize the incoming hydrogen. 
When £^61 is large, 77/ is of order unity, and cosmological effects have begun to govern the growth 
of the blastwave. 
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Fig. 6. — Total ionized mass when ^2 = 1. The horizontal dotted line gives the total 
mass {Mi) ionized by (10^^ erg)^;^! 

in ionizing photons with mean energy (13.6eV)£'Ryd. An 
adiabatic blastwave initially containing (10^^ erg)£^6i ultimately encompasses a mass Masy. The 
diagonal dotted line shows how this asymptotic mass varies with the parameter combination 
-E'6i^5o(l + ■2i)^'^^^iGM ™ ^ universe. The solid line gives the amount of mass (MasyT/f) within 
an adiabatic blastwave when the shock ceases to ionize the incoming gas efficiently. The maximum 
amount of collisionally ionized mass per unit input energy drops with increasing energy because 
high-energy blastwaves transfer most of their energy to the coUisonless component before ionization 
losses become important. Cosmological blastwaves do not ionize the IGM unless their sources emit 
relatively few ionizing photons into intergalactic space. 
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Fig. 7. — Constraints on porosity. COBE limits on tlic y-distortion of the microwave background 
constrain the thermal energy density of the IGM at high redshift. The product of the energy per 
blastwave, Eq, and the comoving number density of blastwave sources, nft(l + z)~^, must therefore 
lie below the line marked y (COBE), if J7 = 1, h^o = 1, and zi = 7. Blastwaves driven into 
the IGM by multiple supernovae introduce metals along with thermal energy. If the metallicities 
of Lya clouds 10~^ times solar) are similar to the mean mctallicity of the IGM, the energy 
density injected by supernovae must be significantly lower than the COBE limit. The line marked 
ZiGM < 10~^ indicates the value of £'o^b(l + that would yield an IGM metallicity of 10~^ 
times solar, assuming each supernova produces 1 Mq of metals and Ojgm = 0.01. The lines labeled 
Q = 1 trace the locus of blastwave parameters that give Q = 1 at z = 4 and z = 2 for blastwaves 
that begin propagating at zi = 7. (See Fig. 8 caption for additional information.) 
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Fig. 8. — Constraints on porosity. Same as Figure 7, except r2iGM = 0.03. The metallicity constraint 
allows Q ^ 1 at z > 2 only if blastwavc sources are much more numerous than galaxies and the 
energy per blastwave is <C 10^'' erg. Quasar blastwaves could be metal-free and need not satisfy the 
metallicity constraint; however, the COBE constraint rules out Q{z > 2) » 1 in quasar blastwaves 
if ^^IGM > 10~^. The line marked QUASARS indicates the comoving number density of quasars 
measured at z = 2 by Boyle (1993), multiplied by a factor of 10 to account for the quasar duty 
cycle (quasar lifetime ~ 10^ yr divided by t 10^ yr at z ~ 5). Blastwaves driven by objects this 
uncommon cannot fill the universe by z > 2, while still satisfying the COBE constraint, regardless 
of their energy output. 
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Fig. 9. — Postshock H column density. This figure shows how the column density of hydrogen 
behind a spherical cosmological shock changes with time in a uniform IGM. The curves indicate 
A^H in blasts beginning at redshifts zi = 4 and 7 with initial energies in the range 10^'^ — 10^^ erg. 
We assume r^icM = 0.01, /150 = 1, and $7 = 1. Note iVn oc ^igm^so""- first, A^h grows. Later, as 
the shock begins to merge with the Hubble fiow, iVn drops. 
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Fig. 10. — Postshock H I column density. This figure gives the neutral postshock column densities 
for the same blastwaves described in Figure 9. The shocks are presumed to be adiabatic (not true 
for Eq ~ 10^^ erg at z < 1.5) and strong enough for the postshock density to be four times the IGM 
density. CoUisional ionization governs postshock ionization balance at early times, indicated by the 
dashed lines. Later, photoionization ( J21 = 1) overtakes collisional ionization. The transition from 
dashed to solid lines occurs where photoionization and collisional ionization are equally important. 
Low energy blastwaves have higher A^hi because their postshock temperatures are lower, allowing 
recombination to act more rapidly. Note A^hi oc ^^igm^so^^- 
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Fig. 11. — Conditions under which shocked clouds cool. Adiabatic shocks that encounter clouds Dd 
times denser than the ambient IGM can drive radiative shocks into these clouds if is sufficiently 
large. The heavy solid lines indicate the shock speeds Vg through the ambient medium below 
which shocked clouds cool in a Hubble time at various redshifts z. Here we show cooling loci for 
CI = 1, h^o = 1, and Oiqm = 0.01. For reference, we also show Vs{z) trajectories for blastwaves 
beginning at zi = 4 with Eq = 10^^, 10^^, and 10^^ erg. Because the postshock cooling time is 
roughly proportional to D^^'^, clouds of modest overdensities can cool much more easily than the 
ambient IGM. Thus, cosmological blastwaves can potentially amplify the neutral column densities 
of preexisting Lya clouds by a large factor. 



